Theory of Type-II Superconductors with Finite London Penetration Depth 
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Previous continuum theory of type-II superconductors of various shapes with and without vortex 
pinning in an applied magnetic field and with transport current, is generalized to account for a 
finite London penetration depth A. This extension is particularly important at low inductions B, 
where the transition to the Meissner state is now described correctly, and for films with thickness 
comparable to or smaller than A. The finite width of the surface layer with screening currents and 
the correct dc and ac responses in various geometries follow naturally from an equation of motion for 
the current density in which the integral kernel now accounts for finite A. New geometries considered 
here are thick and thin strips with applied current, and "washers", i.e. thin film squares with a slot 
and central hole as used for SQUIDs. 

PACS numbers: 74.60. Ec, 74.60.Ge, 74.55.+h 



I. INTRODUCTION 

Thjg statics and dynamics of the Abrikosov vortex 
lineaj and of the two dimensional (2D) pancake vor- 
tices in layered superconductorsB recently was formu- 
lated within continuum approximation in terms of an 
equation of motion of the current density j inside the 
superconductor.0 The resulting integral equation for the 
space and time dependent j(x, y, z, t) contains the as- 
sumed constitutive law for the electric field E„ caused by 
the motion of vortices. In general, the law E = E^Q, B) 
depends on the current density j and on the local induc- 
tion B and describes, e.g., free flux flow, pinning, and 
thermally activated depinning of vortices. 

Besides this, the equation of motion for j depends 
on the geometry of the given problem. Various geome- 
tries have been considered: Thin long strips Jj thin circru 
lar disks and rings Jj thin rectangular platelets or films, 
thick strips]3 and thick circular disks or short cylinders, 
always in a magnetic field applied perpendicular to the 
plane of the superconductor. Further numerical solutions 
for the statics of thin film superconductors of various, 
shapes within the Bean model are given by Prigozhin.El 
In all these geometries, the current density j or the sheet 
current in films J(x, y,t) — (J x ,J y ) — J j(x,y, z,t) dz 
have only one component (flowing along the strip or cir- 
culating in the disk) or may be derived from a scalar po- 
tential g{x,y,t) according to J x — dg/dy, J y — —dg/dx 
since V ■ J = 0. The equation of motion describes thus a 
scalar (j or g) which depends on only one or two spatial 
coordinates. 

A more general continuum formulation including also 
transport currents applied by contacts, applying in prin- 
ciple to arbitrary three dimensional (3D) geometry, and 
accounting also for a finite lower critical field H c \ and 
for a general reversible thermodynamic field H(B), is 
presented in Ref. [lO|, where itj-is applied to the prob- 
lem of the_ geometric barrier .til The elegance of the 
formulationflQ'Ej in terms of an equation of motion for 



the current density is that it explicitly accounts for the 
applied magnetic field or transport current, which in 
other formulations enter via a separate boundary con- 
dition for a differential equation. Our integral equation 
automatically accounts for the inhomogeneous magnetic 
field outside the superconductor, without any need to cal- 
culate this infinitely extended field explicitly during the 
time integration since all spatial integrations extend only 
over the finite specimen cross section. The applied field 
and applied current enter via the screening currents that 
flow on the specimen surface, the information about the 
geometry is contained in a pre-calculated integral ker- 
nel, and the properties of the superconductor enter via 
one single constitutive law E = E v (j,B) if B = /JoH 
may be assumed, or via two laws if some nontrivial law 
B{H) ^ fj, H is used. nnn 

So far, this continuum descriptionBO't2l considered the 
limit where both the vortex spacing a and the London 
penetration depth A are assumed to be much smaller than 
all other relevant lengths. There are, however, two situa- 
tions where a finite London depth A should be considered: 

(i) In bulk superconductors, at small inductions B the 
Meissner state should follow in the limit B — > 0. So 
far, with the above assumption E = E v (j,B) one has 
E = for B — > if the free flux flow law or similar laws 
E v ex B are used. In real type-II superconductors, how- 
ever, one may have E ^ in a surface layer of thickness 
A if j varies with time, since one has E = ^X 2 dj /dt in 
the Meissner-London limit (this London equation means 
that the electric field accelerates the Cooper pairs). 

(ii) In thin films with thickness d, both cases d > A 
and d < A should be described correctly by a continuum 
theory. In particulai_for d <C A the effective penetration 
depth of thin filmsata A — \ 2 /d may become microscop- 
ically large and comparable with the film width. 

In the present paper our previous.-cpj.iti nuum electro- 
dynamics of type-II superconductorsQTJE3 is generalized 
to allow for an arbitrary London penetration depth A. 
The resulting integral equations for j in various geome- 
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tries correctly describe all limiting cases and are easily 
solved numerically. With finite A the stability and speed 
of the numerics is even improved since this A provides a 
well defined inner cutoff length, whereas with A = the 
inner cutoff depends on the chosen grid spacing. 

In the following sections, equations of motion for the 
current density in type-II superconductors with finite 
London depth A are derived for various geometries. For 
illustration, some selected numerical results are presented 
for each geometry, namely, current and magnetic field 
profiles, magnetization curves, ac susceptibilities, and the 
current stream lines in thin films. 



II. VOLTAGE-CURRENT LAW 

In type-II superconductors within continuum approxi- 
mation the local electric field E is composed of two parts, 
originating from the vortex motion and from the Meiss- 
ner surface currents, 



E — E„ 



E 



M ■ 



(1) 



The first term E^jjB) in general is a nonlinear func- 
tion of the current density j due to pinning and thermal 
or quantum depinning (creep). It may be anisotropic 
(a) if the Hall effect of vortex motion is accounted for, 
or (b) if the superconductor is anisotropic, or (c) if j 
is not perpendicular to B. This latter case is still lit- 
tle understood but should occur in most really three 
dimensional geometries, e.g., in a superconductor cube 
with applied magnetic field; one possible model descrip- 
tion in this case uses two critical current densities j c ± 
and j c || for the current components perpendicular and 
parallel to B. Though the parallel current density j|| 
does not exert a Lorentz force on the jirptices it nev- 
ertheless_can trigger a helical instabilitytj'EJ and vortex 
cutting,E£l which leads to a finite j c u. Furthermore, as 
shown by Gurevicht3 some combinations of nonlinearity 
and anisotropy in voltage-current models may lead to 
unstable current distributions. In the following geome- 
try examples, for transparency I shall consider isotropic 
E„ since the extension to anisotropic E„ is straightfor- 
ward. Isotropic voltage-current laws have the form 



(2) 



where p v (j,B) is the resistivity caused by vortex mo- 
tion, and j — |j|, B = |B|, with j 1 B. In the free 
flux flow limit (zero pinning or high current density) one 
has p v — Ps ~ {B / B C 2)p n where B C 2 is the upper crit- 
ical field and p n the normal resistivity just above B C 2- 
For thermally activated depinning many experiments and 
theoretical models yield a logarithmic activation energy 
U(j) = U hx(j c /j), which gives E v = E c exp(~U/kT) = 
E c (j/j c ) n with creep exponent n = Uo/kT ^ 1 and (ar- 
bitrary) "threshold criterion" E c = E(j c ). With Eq. (2) 
this means p v = pAjjjcf where a = n — 1 1. A 



simple model which coratunes this creep model and the 
free flux flow limit readalj 



Pv(j,B) 



Pn_ B U He 
B, 



c2 



i + (J/jcY 



(3) 



The second term in Eq. (1) is caused by the current 
flowing in the Meissner state and is given by the dynamic 
London equation, 



E 



M 



MoA : 



dt 



(4) 



This term describes the acceleration of the massive 
charge carriers (Cooper pairs) by an electric field, 
dj/dt oc E. Thus, a finite London depth A in priucir 
may be introduced into our continuum descriptior 
by using the full voltage-current law, 



E = E„(j,B)+ Mo A^ 



(5) 



This replacement works well in the case of linear re- 
sponse, where it yields a frequency (w) dependent com- 
plex resistivity 



p{u) = p s + iup X 2 



(6) 



in the absence of vortex pinning. In this case the linear 
integral equation for j reduces to an eigen value equation 
from which the complex ac susceptibility of the supercop-, 
ductor for a given geometry follows as an explicit sum.tll 
When E v (j,B) is nonlinear, the integral equation for 
j(jc,t) has to be solved by numerical integration over 
time. It turns out that this time integration becomes 
unstable when the modified voltage-current law (5) is 
inserted. However, when the term containing A 2 is incor- 
porated into the integral kernel as shown in Set. 3, then 
the numerics remains stable and even becomes more sta- 
ble and faster than it was with A = 0. 



III. THICK STRIPS WITH FINITE A 

To fix ideas we first derive the modified integral equa- 
tion for the geometry of an infinite strip or bar with rect- 
angular cross section filling the volume —a < x < a, 
-b < y < b, -L/2 < z < L/2, L > a,b, see Fig. 1. 
In the general 3D geometry, with r = (x,y, z), the 
vector potential Aj(r) caused by the current density 
j(r) = — ^ " 1 V 2 Aj in the gauge V • A = is given by 



A j (r)= A4o /dV 



j(r') 



An r - 



(7) 



with the integral taken over the volume where the cur- 



rent flows. Indeed, using V \v — r' 



-47n5 3 (r - r') 



(<$3 is the 3D delta function) one verifies that Eq. (7) 
yields V 2 Aj(r) = — /i j(r). For a long strip or bar, both 
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j = j(x, y) z and A = A(x, y) z are directed along z. Inte- 
grating Eq. (7) over z and writing from now on r = (x, y) 
we obtain 



A J -(r)= M0 /dVj(r')Qbar(r,r') 



with the 2D integral kernel 

1 L/2 

Qbar(r,r') = — asinh-j - 

2tt r — r 



^— In i r 

2tt r - r' 



(8) 



(9) 



The integral (8) is over the rectangular strip cross sec- 
tion, —a < x < a and — b < y < b, but actually Eq. (8) 
applies to strips with cross sections of any shape. If the 
current distribution is symmetric the integration (8) may 
be restricted to one quarter of the rectangular cross sec- 
tion, e.g., 

Aj(x, y) = Mo fdx' [dy'j(x', y') Q s *™(x, y; x', y 1 ) . (10) 
Jo J o 

For strips with a magnetic field H a applied along y, one 
has j(x, y) = -j(-x, y) = j(x, -y) = -j(-x, -y), thus 



1 in {x t + y2 - ){x i + y i ) (11) 



^bar ^bar ^ - ^_ + y 2_ ) ^2 + y 2 ) 



with x± — x ± x' , y± — y ± y' . For strips with trans- 
port current I a along z one has j(x,y) = j{—x,y) = 
-y) = j{-x, -y), thus 



QZ- = QL 



L 8 



4tt {x 2 _+y±){xt +yV){x\ +y 2 _){x\ +yl) 



(12) 



Note that the strip length L has dropped out in Eq. (11) 
but not in Eq. (12). Therefore, some electrodynamic 
properties of long strips with applied current (e.g. their 
self induction) depend logarithmically on the strip length 
L, while for strips in a magnetic field usually the limit 
L — * oo may be taken. Strips with oblique applied field 
H a = (H ax , H ay ), or with both applied H a and I a , have 
a lower symmetry and the integration (8) then has to be 
taken over the half or full cross section of the strip rather 
than over a quarter. 

To obtain an explicit dynamic equation for the current 
density j{x, y, t) one has to incorporate the applied vector 
potential A , which for the strip may be chosen along z, 
A a = A a (x, y, t) z. In general, A a (x,y,t) may have two 
parts originating from an applied perpendicular magnetic 
field or induction B a (£) = fioH a (t) — (B ax , B ay ) and 
from an applied electric field E a = E a (t) z which drives 
the transport current I a , 



A a (x,y,t)^A^+A^ 



Vx(4fz), E a 



-A 



(13) 



where the dot denotes the time derivative. For example, 
with B a \\y one has A^ = —xB a , and with E a (t = 0) = 



one has = - J Q * E a {t')dt' . Writing the total vector 
potential as A — Aj + A a we obtain from Eq. (8), 

Mo / dVj(r',i)Q b ar(r,r') = A 3 (r,t) =A — A a , (14) 



For the strip the induction law B = — V x E yields 
E = —A. Inserting E from Eq. (5) we have for strips 



A=-E v (j,B)- f i X 2 dj/dt. 



(15) 



thus the r.h.s. of Eq. (14) may be written as 

Aj(r,t) = -J E v (j,B)dt- moA 2 jM) -A a (r,t). (16) 

Inserting this in Eq. (14) and shifting the term MoA 2 j(r, t) 
to the left under the integral we obtain 



Mo J d 2 r' j(r', t)[ Q bar (r, r') + A 2 «5 2 (r - r') ] 

= - JE v (r,t)dt-A a (r,t), (17) 

where 62 (r) = S( x )^(v) 1S the 2D delta function. Tak- 
ing the time derivative, pf, Eq. (17) and introducing the 
inverse integral kerneli3't2l 

K (r, r') = [Q bar (r, r') + A 2 <5 2 (r - r')] _1 . (18) 

we arrive at the equation of motion for j(r,t): 

dj{r,t) 



= -fio 1 d 2 r'K(r,r')[E v (j,B)+A a (r\t)} (19) 



dt 



with A a (r',t) = ~x'B ay (t) + y'B ax (t) - E a (t) and with 
Qbar from Eqs. (9), (11), or (12). In Eq. (19) only the 
electric field caused by the vortex motion E v enters ex- 
plicitly. The London length A enters the inverse kernel 
K(r, r') defined by 



d 2 r"K(r, r")[ Qbar(r", r') + A 2 ^ 2 (r" - r') ] 

= J 2 (r - r') . (20) 



The kernel K(r, r') may be computed by a matrix in- 
version as follows. First, a spatial grid r; = (xi,yi) is 
chosen with appropriate weights wi such that the inte- 
grals over the strip cross section are well approximated 
by a sum, 



Jd 2 rf(r)^Y,^ r ^ 



(21) 



Then we express the definition (20) by such a sum, 

K u ( Qu' + A2 % ) = *y . ( 22 ) 
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where Qij = Q(r,, r^) and Sij equals 1 if i = j and else. 
Solving Eq. (22) for the matrix i<Cy one finds 



Kij = ( Qij Tw i 



+ x 2 s tJ )" 



(23) 



The accuracy of this method is considerably increased by 
choosing a nonequidistant grid with narrow spacing near 
the specimen surface and by taking appropriate diagonal 
terms Q^ ar as described in the Appendix of Ref . 8a. From 
Eq. (23) one sees that finite A 2 increases the (positive) 
diagonal terms of the matrix to be inverted; this makes 
the matrix inversion more stable. 

For numerics we need the equation of motion (19) for 
j(r, t) in discrete form, 



dt 



-Mo 



(24) 



where the vectors ji(t) — j(ri,t), Bi(t) — B(r.i,t), and 
A a i(t) = A a (n,t) = ~XiB ay (t) + y l B ax (t) - E a (t) are 
functions of the time t. The matrix Kfj, Eq. (23), is in- 
dependent of time and has to be computed only once for 
a given geometry and given A. Equation (24) is easily in- 
tegrated over time t starting with ji (t = 0) = and then 
switching on the applied fields B a and/or E a . The result- 
ing magnetic moment per unit length m(t) = (m x ,m y ) 
and total transport current I a (t) (along z) are then ob- 
tained as integrals over the strip cross section, 



m(t) 



dx dy(-xy + yx)j(x,y,t) 



~ + yxi) jj(t) w t , 

i 

/a rb 
dx dyj(x,y,t)rts^2ji(t) 
-a J— b 



IV; 



(25) 
(26) 



Note that the contribution to m(t) of the U-turn of the 
currents at the strip ends (integrals over the x and y com- 
ponents of j, amounting to exactly ^m) is already con- 
sidered in Eq. (25). Note further that, though in exper- 
iments usually the applied current I a (t) is imposed, the 
theory considers a spatially constant electric field E a (t) 
along z which drives the current I a (t) that results from 
the calculation. 

Figures 2-6 show the current density j(x,y) and the 
magnetic field lines of B (i.e. the contour lines of A) for 
strips with aspect ratio b/a = 0.4 in various cases, cf. 
Fig. 1. In Fig. 2 the strip is in the Meissner state, i.e., no 
vortices have penetrated. Shown is the screening current 
density for finite London depth A = 0.025a for a strip in 
applied field (top) and for a strip with applied current 
(bottom). Similar figures for strips with quadratic cross 
section are depicted in Ref. [t9| Note the sharp peak of 
j(x,y) in the four corners, which has finite height. This 
high local current density favors the nucleation of vortex 
quarter loops from the corners of the strip when the ap- 
plied magnetic field or current exceed a certain threshold. 



Figures 3 and 4 show current density and field lines 
for strips in an increasing applied field H a for two values 
of the London depth A/a = 0.025 and 0.1. Since the as- 
sumed voltage-current law is very steep, E v oc j n with 
n = 101, a flat saturation of j(x, y) occurs at the critical 
value j c , like in the Bean model. The field of full penetra- 
tion of this strip ia3 H p = (j c b/n)[(2a/b) arctan(6/a) + 
ln(l + a 2 /b 2 )] = 0.4945aj c for b/a = 0.4 and H p = 
2.52 • (2bj c ) for b/a = 0.001. 

Figures 5 and 6 show the same strips but with increas- 
ing applied current I a and with no field applied. In this 
case, when the critical current I c — ^abj c of the strip is 
reached, vortex rings penetrate continuously and annihi- 
late in the center of the strip. With further increasing 
I a > I c , the dissipation and voltage drop increase steeply. 
More results for strips with both applied field and current 
will be published elsewhere. 

Figures 7 and 8 show magnetization loops m{H a ) of 
strips in perpendicular applied field H a (t) = Hq sin Lut at 
two amplitudes Hq, for five values of the London depth 
A, and for two creep exponents n entering in E v (j) — 
E c (j/j c ) n . Figure 7 is for a thick strip (b/a = 0.4) and 
Fig. 8 for a thin strip (b/a = 0.001). With our dimcn- 
sionless units a — j c = E c = (j,q = 1, the used circu- 
lar frequency oj = 1 corresponds to w = E c /(^j, Q j c a 2 ) in 
physical units. Note that with increasing A the hysteresis 
loop becomes more narrow and finally collapses to a line, 
i.e., the magnetic response becomes reversible. This is 
so since with increasing A the screening current density 
decreases and can no longer depin vortices, except when 
H a is large or n is small. 



IV. THIN STRIPS 

This section considers the limit of thin strips, b -C a, 
with B a applied perpendicular (along y) and I a applied 
parallel (along z) to the strip. In this limit only the cur- 
rent density integrated over the film thickness matters, 
called sheet current and directed along z (like j), 



J(x,t)= / j(x,y,t)dy. (27) 
J-b 

Integrating Eq. (17) over y we obtain (see also Ref. 
jiQ j dx J(x , t) 



l ln ]x^\ +M{x - X,) 



- / E v (x,t)dt - A a (x,t) 



(28) 



where A = A 2 /d is the effective penetration depth of 
thin films with thickness d — 2b < X. Here I have used 
the approximate constancy of the current density along 
y, j(x,y) « J(x)/d. Initially it was not clear to mc if 
thin film expressions of the type (28) describe also the 
dynamics of superconductor strips and not only the stat- 
ics, which was successfully considered, e.g. in the static 
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Bean model .calculations of thin disks ,0 thin-strips ,EHl~E£l 
and ellipses rJ But recently we have provenO that not 
only j(r, t) but also the electric field E(r, t) is practically 
constant over the film thickness; deviations from this con- 
stancy occur only near the (penetrating or exiting) flux 
front but are restricted to a transverse length scale of or- 
der d. This result was obtained in the limit A — ► 0, but 
it applies all the more for finite A. Thus, the dynamics 
of the sheet current J * jd is well described by Eq. (28). 

Inverting Eq. (28) and taking the time derivative we 
obtain the equation of motion for the sheet current in 
thin strips, 



(Jl j(x,t) 



dx'K(x, x')[E v (x', t) + A a (x', t) ] (29) 



with the inverse integral kernel 



K(x,x') = 



2tt x 



L 



AS(x - x') 



(30) 



and with E v (x' ,t) = E v {j,B) depending on j — 
J(x',t)/d and B = \B(x',t)\, and with A a (x',t) = 
—x'B a (t) — E a (t). For strips with either applied field 
B a or applied current I a , the integration in Eq. (29) may 
be restricted to half the strip width, < x < a. For 
B a ^= 0, I a = one has J{x) = —J{~x) and a symmetric 
kernel [cf. Eq. (11)] 



K B (x,x') = 



2tt \x 



A6(a 



(31) 



For B a = 0, I a y£ one has J(x) — J(—x) and a sym- 
metric kernel [cf. Eq. (12)] 



Kj(x,x') 



1 



In 



L' 2 



2tt \x-x'\(x+x') 



-+AS(x-x') 



(32) 



Note that the strip length L has dropped out in the kernel 
(31) but not in the kernel (32). The complex resistivity 
of thin strips thus in general depends on the logarithm 
of the strip length. 

Figures 9 and 10 show the profiles of the sheet current 
J(x) and perpendicular induction B y {x) of thin strips 
with increasing applied magnetic field (Fig. 9) or current 
(Fig. 10) for various values of the effective penetration 
depth A = \ 2 /d; d = 2b is the strip thickness and 2a 
the strip width, here b/a — 0.001. Note that for larger A 
the profiles J(x) become smoother, but the bend of J(x) 
at the point where J starts to deviate from the critical 
sheet current J c = dj c remains sharp. For A/a > 0.2, the 
sections of J{x) where J\ < J c are almost straight lines 
in Fig. 9 and almost parabolas in Fig. 10. 

Figure 11 shows J(x) and B y (x) in thin strips which 
are exposed to a high magnetic field H a ^> H p before 
a current I a is applied, ranging from to the criti- 
cal current I c — 4a6j c . First, when I a = 0, one has 
J(x) w J c sign(x), and finally, when I a rj I Cl one has 
J(x) w J c (equal signs would apply in the Bean limit 



n — > oo). At intermediate currents I a < I c , in the half 
strip < x < a one has J w J c = constant, and in 
the other half —a < x < the profile J(x) is symmetric 
about the point x = —a/2. For A = 0, n>l, the known 
Bean profiles for thin-strips with transport current apply 
to this half strip £313 



V. AXIAL SYMMETRY 

The above method for incorporating finite A can also 
be applied to axially symmetric problems, i.e. to super- 
conductors with an axis of rotational symmetry in a mag- 
netic field applied parallel to this axis (along y). A sim- 
ple example are disks or short cylinders with radius a and 
thickness d = 2b, but also other shapes like rings, toruses, 
cones, ellipsoids, etc. are easily described by introducing 
a y dependent radius a(y). In all these cases the current 
density and vector potential have only one component 
directed along the azimuthal unit vector <p, j = j(r,y)(p 
and A = A(r,y)(p, with r = (x 2 + z 2 ) 1 / 2 . Integrating 
the 3D Eq. (7) over the angle tp — arctan(z/a;) we obtain 
the vector potential Aj caused by the current density j 
circulating in an axially symmetric conductor, 



r b l-a(y') 

A J (r,y) = f i I dy' / dr'j(r',y') 



-b 



t(r,y,r',y') (33) 



with the kernel! 



f(r, r',rf) 



Q^(r,y;r',y') = f(r,r',y-y'), (34) 
^dip —r' cos ip 



o 2tt (rj 2 + r 2 + r' 2 — 2rr' cost/?) 1 / 2 



(35) 



When the conductor has a symmetry plane at y = 0, 
then the integration over y' in Eq. (33) may be restricted 
to < y < b if the kernel Q ax is replaced by a symmetric 
kernel, 

Q cyl (r, y; /, y') = f(r, /, y - y') + /(r, r', y + y'). (36) 

If in addition the radius is constant (like in disks and 
cylinders), Eq. (33) simplifies to 



Aj{r, y) = Mo / dy' / dr'j(r', y') Q cyl (r, y; r> , y 1 ) . (37) 
Jo Jo 

Comparing Eqs. (33) and (37) with Eqs. (8) and (10) 
for the strip or bar, one easily verifies that the analogue 
to Eq. (17) for axial symmetry is identical to Eq. (17) 
but with Qbar(r,r') replaced by Qax(r,r'), Eq. (34) or 
Qcyi(r, r'), Eq. (36), where now r = (r, y), r' = (/, y'). In 
axial symmetric problems no transport current and thus 
no electric field is applied but only a magnetic field or in- 
duction B a = V x (A a (r, y, t)<p) = (B ar , B ay ). In general 
this applied field may be inhomogeneous, e.g. when levi- 
tation forces are to be computed. If B a is homogeneous 
one has A a {r,y,t) — —%B a (t). With this notation, for 
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the general axially symmetric geometry the equation of 
motion for the circulating current density is identical to 
Eq. (19) but with a different inverse kernel 



^(r,r') = [Q a ,(r,r')+A 2 5 2 (r-r')]- 1 , 



(38) 



with Q ax from Eq. (34). If the superconductor has a 
symmetry plane at y — and the applied field is homo- 
geneous (or symmetric about the plane y = 0) one may 
replace Q ax in Eq. (38) by Q cy \ from Eq. (36), 



K(v,v') = [Q cyl {v,v')+\H 2 (v-v')Y 



(39) 



These inverse kernels K (r, r') may be computed by the 
same method discussed below Eq. (20) if the correct 
integration area (or spatial grid) is used as defined in 
Eqs. (33) and (37). For the best choice of the diagonal 
terms Qa see Appendix of Ref. 8a. 

The axial magnetic moment of a disk or any other ro- 
tationally symmetric conductor is 



m(t) 



dy 



a(y) 



dr r 2 



j{r,y) . 



(40) 



If the applied field is periodic, H a (t) = Ho sin Lot, the 
(in general nonlinear) complex susceptibilities ofrthe disk 
may be defined as \v = Xu ~ *x", v — 1,2,3.. 



irH, 



o Jo 



m(t) e" ll/UJt d{ujt) . 



(41) 



Here \i is the fundamental susceptibility and the Xv with 
v > 1 correspond to higher harmonics, which are absent 
for linear response. The Xv usually are normalized such 
that for Hq — > or uj — > oo the ideal diamagnetic sus- 
ceptibility xi(0, w) = —1 results. This is achieved by 
dividing all x»i Eq. (41), by the absolute value of the 
initial slope [dm(H a ) / 'dH a ]n a =o- 

As one example, Fig. 12 shows the real and imagi- 
nary parts of the nonlinear fundamental susceptibility 
Xi{H q ,uj) — x — x" ~ °f a thick disk with b/a = 0.5 
plotted versus the ac amplitude Hq for constant fre- 
quency u> = E c /(poj c a 2 ), creep exponent n = 11, and 
for various London depths A/a = 0.025 ... 1. The same 
data are depicted in Fig. 13 as polar plot, x" versus — 
In both presentations x(^o) sensitively depends on the 
parameters n and A/a and on the geometry. More data 
Xv(Hq,u>) are available from the author. 



VI. THIN PLATES AND FILMS 

The geometry of thin plates or planar films with arbi- 
trary shape in a perpendicular magnetic field i? ||z differs 
from the geometries of Sets. 3 to 5 in that now the cur- 
rent density is no longer a scalar but has two components, 
j(x,y, t) = (j x ,jy)- However, because of the strict rela- 
tion diyj = 0, this planar j may be derived from a scalar 



potential (or magnetization, stream function) g(x,y,t). 
Since we are interested in the thin film limit we consider 
the sheet current 3(x,y,t) = (J x ,J y ) — J j(x,y, z,t) dz, 
Jx = dg/dy, J y = —dg/dx. The stream lines of 3(x,y) 
coincide with the contour lines g(x, y) = const. Thus, one 
may choose g(x, y) = on the edge of the film since the 
current flows along the edge. In this section I generalize 
the theory of Ref. g to finite London depth A. 

The magnetic field H = H z (x,y)z in the film plane 
2 = is related to the local magnetization g(x, y) by an 
integral over the film area 5, 



H z (r) = H a + fd 2 r' Q film (r, r') g(r') 
Js 

with r = [x, y) and r' = (x 1 ', y'). The integral kernel 

2z 2 - p 2 



Qfiim(r,r') = lim 



2^0 47r(z 2 + p 2 ) 5 / 2 



(42) 



(43) 



with p 2 — (x — x') 2 + (y — y') 2 , gives the magnetic field 
generated by a magnetic dipole of unit strength posi- 
tioned in the plane z = at (x',y'). Explicit expressions 
for films with rectangular shape are given by Eqs. (42) to 
(46) of Ref. ^ in form of Fourier series. To incorporate a 
finite London depth A we write the voltage-current law 
of the film in the form of Eq. (5) but now as a function 
of the sheet current J(x,y) = j(x, y)d(x, y). We consider 
first a uniform isotropic film with constant thickness d. 
The local electric field is then 



E(J,S) =p s (J,5)J(r,i)+ Mo Aj(r,i) 



(44) 



where p s = pjd is the sheet resistivity and A = A 2 /d 
is the effective magnetic penetration depth of the film. 
This electric field is related to the induction B = poll 
by the induction law B = —V x E, which in the film plane 
z = means B z = dE x /dy-dE y /dx. With Eq. (44) and 
Jx — dg/dy, J y = —dg/dx, we obtain 

B z (r, t) = y[ Ps Vg(r, t)] + /i AV 2 .g(r, t) . (45) 

Taking the time derivative of Eq. (42), inserting B z = 
PoHz from Eq. (45), and combining the two terms con- 
taining g, one arrives at 

/ dV [Q fi im(r, r') - 5 2 (r - r')AV 2 ] g(r' , t) 
Js 

= f(r,t)-H a (r,t), (46) 



with 



/(r,t)= M o 1 V[p s (r,t)Vp(r,t)] 



(47) 



Inverting this we obtain the equation of motion for 

g(x,y,t): 



g(r,t)= d 2 r'K{v,v')[f(v',t)-H a (v',t)] (48) 
Js 
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with /(r, t) from Eq. (47) and with the inverse kernel 

K(r,r')= [Q mm (r,r')-S 2 (v-r')AV 2 }- 1 . (49) 

This kernel K{ri,Vj) may be evaluated on a grid as a 
Fourier series with discrete k vectors K. If the Fourier 
coefficients of Qahn( r *> r j) are Qkk', given by Eq. (46) 
of Ref. (6) for a rectangular platelet, then the Fourier 
coefficients of K(ri, rj) are the inverse matrix 



[Qkk'+A^kk'T 1 . 



(50) 



Figure 14 shows the current stream lines in a thin rect- 
angle with side ratio b/a — 0.7 for two effective penetra- 
tion depths A/a = and A/a = 0.1 and three applied 
fields H a /J c = 0.001, 0.5, and 1.5 (ideal screening, par- 
tial penetration, full penetration). Note that finite A 
rounds the sharp corners of the rectangular stream lines 
and slightly delays vortex penetration. 

If the superconductor film is nonuniform [e.g. has spa- 
tially varying thickness d{x,y)\ and/or anisotropic [e.g. 
has two (in general nonlinear) resistivities p xx and p yy 
and/or two London depths X x and X y ], then a rather 
general current- voltage law E(J,£?,r) = (E x ,E y ) is 

E X = P SX J X + P()A X J X , Ey — PsyJy + P()AyJy , (5l) 

with the two sheet resistivities 

Pxx (3,B,v) _ p V y(3,B,r) 

Ps.\v)- d{r) , Ps V (r)- ^ , (52) 

and two effective penetration depths 



A(r) = MW A (r) = W 
x[ ) ~ d(r) ' y( ' ~ d(r) ■ 

The generalized Eqs. (45) and (46) are then 



z ,\ d ( dq . dg 

d ( dg dg 
- — n„ - — 



(53) 



Of og og \ 



d 2 r l Q(r,r')g(r , ,t) = f(r : t)-H a (r,t), 



(54) 



(55) 



with 



Q(r,r') = Qfiim(r,r')-£ 2 (r-r') 

( ® A ^ ^ A ^ ^ 

\dx y dx dy x dy) ' 

oi \ -i T d ( dg\ d ( dg\ 

= ^ (/'-'< ' IhjV'-Ui,) 



(56) 
(57) 



The equation of motion for g{x 1 y, i) now is still given by 
Eq. (48) but with /(r,t) from Eq. (57) and with the in- 
verse kernel K(r,r') — Q(r,r') _1 , Eq. (56). The kernels 
Q and K may be evaluated by Fourier transformation. 



Some figures of current stream lines and field profiles for 
thin rectangles with anisotropic pinning and A = are 
depicted in Ref. [2^. 

When A = 0, tho-kiverse kernel K = Qg^ can be ob- 
tained by iterationjH3 see also the Appendix of Ref. |27|: 

'- ^(r-v') + 

IS 



K(r, r') 



'C(r) 



K(r",r')-K(r,r>) d2r „ 



47T r - 



r-"|3 



(58) 



starting with K(r,r') = 0. Here C(r) is defined as an 
integral over the infinite area S outside the film or as a 
contour integral along the film edge: 



2tt 



R{4>) 



(59) 



with R = |r — r'|. For rectangular films < a, \y\ < b) 
one has explicitly 



c ( x > v) = J2 [( a _ p x ^ 2 + ( 6 - w)" 



(60) 



with p,q = ±1. Expression (60) may also be used for 
rectangular films with holes or slots if the integral in 
Eq. (58) is restricted to the area of the real film. 

For numerics one needs the discretized versions of 
Eq. (42) and its inverse. Introducing a 2D grid r^ = 
(xi,yi) with weights to, [J2i w i — °f- Eq. (21)] and 



writing Q l} 
H a , and gt = 



H z {vi) 



- g(vi) one has 



(61) 



The inverse kernel Kij = (QijWj)~~ l cannot be obtained 
directly by discretizing Qfii m (r,r'), Eq. (43), but it may 
be obtained by iterating Eq. (58). This iteration can be 
solved in one step by inverting a matrix, 



Sij \Ci + ^ w lQil ) - (1 - Si^WjQij 



(62) 



where (7j = C(rj) and = l/(47r|rj — rj\ 3 ). Note 
that the terms in Eq. (62) must not be rearranged since 
qu = oo. 

Using this matrix K ^ the current stream lines in ide- 
ally screening flat films of arbitrary shape or of films with 
trapped vortices are easily calculated. If the applied field 
is constant, hi = —H a , the stream lines are the contour 
lines of the trace of matrix Ka: 



ij ■ 



(63) 



If a vortex with flux $o is trapped at a grid point rj 
in the film, the magnetic field is H z (v) = $>oS(r — rj) s» 
(&o/wj)5ij at r = rj, since we assume here A = 0. Thus, 
gi is just one row of the matrix: 



9i = ($o/wj) ■ K v 



(64) 
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In general, the stream function gi is the linear superposi- 
tion of the contributions of the applied field and trapped 
vortices, and possibly of currents applied by contacts. 

The following relations are useful: The difference gi—gj 
equals the total current crossing any line connecting the 
points Yi and Yj; the function ±.g(r) is the potential in 
which a probe vortex (or flux bundle) of appropriate sign 
moves; the Lorentz force on a vortex is perpendicular to 
the contour lines of g{r); and the matrix Kij is propor- 
tional to the interaction energy between two vortices at 
Yi and Yj. All these statements apply to any film shape 
and to any A if the general kernel K, Eq. (49), is used. 
For A = 0, the kernel K is given by Eqs. (58) or (62). 

Thus, a vortex will nucleate at the position on the film 
edge where the stream lines are densest, and then move 
to an cxtremum of g(r) if it is not pinned by material 
inhomogeneities. 

Figure 15 shows the so called washer geometry that 
is used for SQjUIDs (Superconducting Quantum Interfer- 
ence Devices) ,E3 namely, a superconductor thin square or 
rectangle (or circular disk) with a slot and central hole. 
See Refs. |3ij-|34| for theories of such washers with finite 
effective penetration depth A. After nucleation at one 
border of the slot, a vortex without pinning will move 
towards a maximum of g(r). For the square washer in 
Fig. 15 (top left) there is one flat maximum of g(r) to 
the left of the central hole. For the rectangular washer 
(bottom) there are three such maxima, to the left of the 
center and above and below the slot, which are sepa- 
rated by two saddle points above and below the central 
hole. A vortex may thus enter near the middle of the 
slot, jump to the nearby maximum of g, and then jump 
over the saddle point to the wide maximum at the left. 
The field-driven thermally activated motion of vortices 
generates low- frequency noise in the SQUID, which may 
be suppressed by introducing pinning centers, e.g. tiny 
holes H 

In the examples of Fig. 15, A = (or A <C a) was as- 
sumed. In this case the vortices interact with the screen- 
ing currents of the applied field and with other vortices 
only via the stray field outside the film. In general, at 
distances r <C A the point vortices in thin films interact 
via a logarithmic potential similar to the interaction of 
line vortices in the bulk. At distances r 3> A the interac- 
tion is mediated by the stray field, is of long range, and 
depends on the shape of the film. This long range part 
of the vortex interaction (or flux-bundle interaction) is 
given by the integral kernel K, Eq. (58) and (62). For 
infinitely extended films the A dependent vortex inter- 
action and the magnetic field m£ a vortex in thin films 
(d <C A) were derived by Pearl Ji3 see also Refs. | 36. For 
infinite films of arbitrary thickness the vortex interaction 
and magnetic field are obtained in Refs. [57 38 



VII. SUMMARY 

In this paper finite London penetration depth A is in- 
troduced into known continuum methods which compute 
the electromagnetic response of type-II superconductors 
in various geometries. In addition, solution methods for 
some new geometries are presented, namely, thick and 
thin strips with transport current and with both applied 
current and applied magnetic field, and thin films of arbi- 
trary shape, e.g., the washer geometry used for SQUIDs. 
The inverse integral kernel K{r, r') or the matrix 
defined in Eqs. (58) and (62) have the simple physical in- 
terpretation of the interaction energy between two point 
vortices in a film of given shape. The stream function 
g(r) of Eqs. (63) and (64) is just the potential in which 
a probe vortex (or flux bundle) in the film moves and 
which is caused by the applied perpendicular magnetic 
field (and/or applied current) and by the other vortices. 
In the depicted case of small A = A 2 /d, this interaction 
is mediated only via the magnetic stray field outside the 
thin film, which crucially depends on the shape of the 
film. For finite films with arbitrary A, the vortex inter- 
action is implicitly given by the inverse kernel K(r,r'), 
Eq. (49). Explicit A dependent expressions for concrete 
film shapes will be given elsewhere. 
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FIG. 2. A strip with aspect ratio b/a = 0.4, cf. Fig. 1, in the 
Meissner state with London penetration depth A = 0.025a. 
Shown is the current density in a quarter of the cross section 
and the magnetic field lines (inset). Top: In perpendicular 
applied magnetic field H a - Bottom: With applied current. 




FIG. 1. Geometry of long superconductor strips with rect- 
angular cross section 2a x 2b in a perpendicular applied 
magnetic field H a (left) or with applied electric current 
I a = I = J j(x, y) dx dy (right). 
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FIG. 3. Current density j(x,y) (left) and magnetic field 
lines (right) of a superconductor strip with aspect ratio 
b/a = 0.4 in increasing magnetic field H a = 0.05, 0.15, and 
0.35 in units aj c (j c = critical current density). The dashed 
lines are contours of the current density at j/j c = ±0.75, 
±0.45, and ±0.15. The superconductor is characterized by a 
pinning-caused voltage-current law E oc (j/jc) 101 and by a 
small London penetration depth A/a = 0.025. 




FIG. 4. As Fig. 3 but for larger London depth A/a = 0.1. 
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FIG. 5. Current density j(x,y) (left) and magnetic field FIG. 6. As Fig. 5 but for larger London depth A/a = 0.1. 
lines (right) of a superconductor strip with aspect ratio 
b/a — 0.4 with applied current I a = 0.22, 0.66, and 0.94 
in units of the critical current I c — 4abj c - The dashed lines 
are contours of the current density at j/j c = 0.75, 0.45, and 
0.15. As in Figs. 3, 4 a voltage-current law E oc (j/jc) 101 was 
used and a small London penetration depth A/a = 0.025. 
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FIG. 7. Magnetization curves m(H a ) of a superconductor 
thick strip with aspect ratio b/a = 0.4 in a perpendicular field 
H a (t) = Hosinuit at two amplitudes Ho/H p = 2 (left) and 
H /H p = 1 (right) with H v = 0.4945aj c the field of full pen- 
etration. Fields in units aj c and magnetic moment m in units 
a 3 j c (per unit length) . Shown are virgin curves and hysteresis 
loops for two creep exponents n = 101 (top) and n = 7 (bot- 
tom) and for London depths A/a = 0.02, 0.1, 0.2, 0.4, and 
0.6. 




-3 -2 -1 ,q 1 23 -1 q 1 



FIG. 8. Magnetization curves m(H a ) of a superconductor 
thin strip with aspect ratio b/a = 0.001 in a perpendicular 
field H a (t) = Ho smut at two amplitudes Ho = 3 (left) and 
Ho — 1 (right) in units of the critical sheet current J c = dj c 
(d — 2b). Magnetic moment m in units a 2 J c (per unit length), 
Hp — 2.52 J c . Virgin curves and hysteresis loops for two creep 
exponents n = 101 (top) and n = 11 (bottom) and five effec- 
tive penetration depths A = \ 2 /d: A/a = 0.02, 0.1, 0.2, 0.4, 
and 0.8. 
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FIG. 9. Profiles of sheet current J(x) (dashed lines) and 
perpendicular induction B y (x) (solid lines) of a thin strip in 
increasing perpendicular magnetic field H a = 0.10, 0.22, 0.35, 
0.49, 0.65, 0.82, 1.0, 1.23, 1.46, and 1.72 (curves from bottom 
to top). J and H a are in units J c — dj c (d = 2b <C a), B 
in units fj,oJ c . For six effective penetration depths A = A 2 /d: 
A/a = 0, 0.05, 0.1, 0.2, 0.4, and 1. Creep exponent n = 51. 



FIG. 10. As Fig. 9 but for strips with increasing applied 
current I a /I c = 0, 0.2, 0.4, 0.6, 0.8, 0.9, 0.965, and 1. 
Ic — 4abj c is the critical current of the strip. 
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FIG. 11. As Fig. 9 but for a thin strip in large magnetic 
field H a ^> Hp to which an increasing current is applied, 
Ia/Ic = 0, 0.2, 0.4, 0.6, 0.8, 0.9, 0.97, and 1 (from bottom to 
top) with I c = 4abj c . A/a = 0, 0.05, and 0.2. n = 51. 
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FIG. 12. The nonlinear complex ac susceptibility 
x(Ho,ui) = x' — i-x" of a type-II superconductor thick disk 
with radius a and half thickness b = 0.5a plotted versus the 
amplitude Ho of the ac magnetic field in units of the Bean 
field of full penetration H v = 0.722aj c . The curves are for dif- 
ferent London depths A/a = 0.025, 0.05, 0.1, 0.15, 0.2, 0.35, 
0.5, 0.7, and 1. E v oc j n , n = 11, w = E c /(noj c a 2 ). 
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FIG. 13. The nonlinear susceptibility of Fig. 12 plotted as 
x" versus — x' ■ 
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A/a = " A/a = 0.1 



FIG. 14. The current stream lines in a thin superconduc- 
tor rectangle (b/a = 0.7) with A/a = (left) and A/a = 0.1 
(right) at applied perpendicular fields (from top to bottom) 
H a /J c = 0.001, 0.5, and 1.5. 




FIG. 15. Thin superconductor square (a = b) and rectan- 
gle (a/b = 3) with slit and hole (washer) in the Meissner state 
with A = 0. Shown are the current stream lines for constant 
applied field (top left, bottom) and for zero applied field but 
with one or two vortices trapped at various positions. 
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